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A general class of nonadditive sticky-hard-sphere binary mixtures, where small and large spheres represent 
the solvent and the solute, respectively, is introduced. The solute-solute and solvent-solvent interactions are 
of hard-sphere type, while the solute-solvent interactions are of sticky-hard-sphere type with tunable degrees 
of size nonadditivity and stickiness. Two particular and complementary limits are studied using analytical 
and semi-analytical tools. The first case is characterized by zero nonadditivity, lending itself to a Percus- 
Yevick approximate solution from which the impact of stickiness on the spinodal curves and on the effective 
solute-solute potential is analyzed. In the opposite nonadditive case, the solvent-solvent diameter is zero and 
the model can then be reckoned as an extension of the well-known Asakura-Oosawa model with additional 
sticky solute-solvent interaction. This latter model has the property that its exact effective one-component 
problem involves only solute-solute pair potentials for size ratios such that a solvent particle fits inside the 
interstitial region of three touching solutes. In particular, we explicitly identify the three competing physical 
mechanisms (depletion, pulling, and bridging) giving rise to the effective interaction. Some remarks on the 
phase diagram of these two complementary models are also addressed through the use of the Noro-Frenkel 
criterion and a first-order perturbation analysis. Our findings suggest reentrance of the fluid-fluid instability 
as solvent density (in the first model) or adhesion (in the second model) is varied. Some perspectives in terms 
of the interpretation of recent experimental studies of microgels adsorbed onto large polystyrene particles are 
discussed. 

I. INTRODUCTION 

Many years ago, Asakura and Oosawa^ (AO) pro¬ 
vided an explanation of the clustering and gelation phe¬ 
nomenon occurring when small nonadsorbing polymers, 
such as polystyrene (PS), were added to a solution 
of large spherical colloids, say polymethylmethacrylate 
(PMMA). The basic idea is illustrated in Fig. 1 consid¬ 
ering two PMMA colloids, modeled as big spheres, im¬ 
mersed in a fluid formed by a uniform background (that 
we will neglect henceforth) as well as by PS particles, 
assumed to be small noninteracting spheres that, how¬ 
ever, experience a hard-sphere (HS) interaction with the 
larger ones. Under these conditions, when the separation 
between the two large spheres is less than the diameter 
of the small spheres (see Fig. 1), there is an unbalanced 
pressure of the “sea” of small spheres, providing an en- 
tropic gain compared to the case when the separation is 
large, that can be reckoned as an effective attractive in¬ 
teraction driving the clustering of large colloidal spheres. 

In real systems, however, the solvent particles do not 
always behave as an ideal gas or interact only sterically. 

Typically, they experience an additional short-range at¬ 
traction (or repulsion) with the solute, usually due to dis- 
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FIG. 1. Cartoon of the AO depletion interaction. The shaded 
region around each solute represents the volume excluded to 
the centers of the solvent particles.^ 


persion forces.The simplest way of accounting for a 
short-range solute-solvent attraction is by means of Bax¬ 
ter’s sticky-hard-sphere (SHS) modeF characterized by 
a stickiness parameter Tsi- Both issues (solvent-solvent 
repulsion and solute-solvent short-range attraction) were 
recently addressed by two experimental studies®’® on ad¬ 
sorbing microgels (MG) to large PS latex suspension. In 
this case, the expected mechanism will be clearly differ¬ 
ent, as illustrated by Fig. 2, inspired by a similar figure 
of Ref. 8. 

Let ai and CTs be the diameters of the large and small 
spheres, respectively, and suppose we fix the volume frac¬ 
tion r]i of the large colloidal spheres and gradually in¬ 
crease the volume fraction rjs of the small solvent spheres. 
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FIG. 2. Different mechanisms occurring in the presence of 
a short-range attraction between solvent and solute, as the 
solvent concentration increases: (a) rjs = 0, HS behavior; (b) 
0 < 7?s < Vs j small fraction of solvent particles act as 
bridges connecting the solute into a cluster; (c) rjs ~ Vs, most 
of the solute colloids are covered and again behave as HSs 
with an effective diameter as + ar, (d) rjs > Vs, the “dressed” 
solutes feel an effective depletion attraction. The dark and 
light shaded regions around the solute particles in panel (d) 
represent the effective solute size and the effective volume 
excluded to the centers of the solvent particles, respectively. 


In the absence of solvent particles, the solute particles 
will behave essentially as HSs, as depicted in Fig. 2(a). 
Now imagine we gradually add the small solvent parti¬ 
cles. Because of the solute-solvent attraction, they will 
tend to get adsorbed on the surface of the larger particles 
and mediate an effective attraction between them. This 
bridging mechanism destabilizes the solution as the large 
colloidal spheres tend to form aggregates, as schemat¬ 
ically illustrated in Fig. 2(b). The global effect is the 
formation of a gel phase caused by a free-energy driven 
phase separation of the large and small spheres. As Vs in¬ 
creases, solvent particles tend to progressively cover the 
solute surface, as depicted in Fig. 2(c). We can easily 
estimate^® the critical value ry* at which all large spheres 
will be completely covered to be 77* ~ ry; (27r/-\/3) cts /cri, 
as discussed in Appendix A. At this point, all the solute 
colloids can be “fully covered” by solvent particles and 
they will behave essentially again as HSs with an effective 
diameter ai + ag, with a few additional free solvent parti¬ 
cles. This situation is pictured in Fig. 2(c). Upon adding 
further solvent particles, however, depletion forces be¬ 
tween the small and the covered colloids set in [Fig. 2(d)] 
and phase separation occurs again, this time entropically 
rather than free-energetically, as in the case of Fig. 2(b). 
A useful way to represent the phase diagram of such a 
binary mixture is through an {vhVs) diagram at fixed 
values of size ratio q = Gsjoi and stickiness In this 
diagram, there will be geometrically inaccessible regions, 
for example for Vs or r]i larger than TrlS-s/^, and lines 
separating the various phase coexistence regions. The 
topology of the phase diagram would be controlled by q, 
while Vs would play the role of an inverse temperature. 

Motivated by these new experimental perspectives, re¬ 
cently Chen et al.^^ considered a HS-SHS binary mixture 
where one can tune the attraction parameter Tsi between 


the unlike spheres, with like spheres only interacting via 
HS interactions. Note that this is the same model already 
studied by Fantoni et al.,^^ as well as by other groups.^® 
The study of Ref. 11 provided a well defined framework 
to rationalize the experimental results obtained in Refs. 
8 and 9. 

In the present work, we will build upon this idea and 
go further to introduce also an additional —and, to the 
best of our knowledge, new— model that has the inter¬ 
esting feature of including the standard AO modeH’^^ as 
a particular case. In both cases, we will illustrate how 
an effective one-component solute-solute interaction po¬ 
tential can be obtained and the merits and drawbacks of 
this procedure. 

Both models can be seen as extreme realizations 
of a general class of nonadditive sticky-hard-sphere 
(NASHS) binary mixtures where the small-small (or 
solvent-solvent) and large-large (or solute-solute) inter¬ 
actions are of HS type with diameters ass and an = ai, 
respectively, while the small-large (or solvent-solute) in¬ 
teractions are of SHS type characterized by a cross di¬ 
ameter a si = {as + ai)/2 = ai{l + q)/2 and an inverse 
stickiness Tsi. Note that here we denote by as= qai the 
diameter of the small spheres as seen by the large ones, 
while ass is the diameter of the small spheres as seen 
by themselves. Thus, the nonadditivity of the unlike in¬ 
teractions is monitored by the ratio ass/as < 1 (where 
we have restricted ourselves to zero or positive nonad¬ 
ditivity). The NASHS class reduces to the nonadditive 
hard-sphere (NAHS) class if the solute-solvent stickiness 
is switched off. 

In the first model that we will study one has ass /= 
1, so that the HS interactions are additive. This model, 
denoted henceforth as the additive sticky-hard-sphere 
(ASHS) model, is the one depicted in Fig. 2 and con¬ 
sidered in Refs. 11 and 12. Interestingly, the ASHS 
model can be solved exactly within the Percus-Yevick 
(PY) approximation^^^^® and the instability region in the 
{Vs,Vi) plane enclosed by the spinodal line can be com¬ 
puted. This will be found to form a closed region, in 
agreement with previous results. 

The second model represents an extreme case of pos¬ 
itive nonadditivity, namely CTss/o’s = 0, i.e., the solvent 
spheres behave among themselves as an ideal gas. This 
particular case of the general class of NASHS models re¬ 
duces to the conventional AO model if the stickiness is 
switched off (i.e., Tgi —t 00 ). Because of that, we will 
term this model as the sticky Asakura-Oosawa (SAO) 
model. The need of supplementing the AO model with 
a short-range solute-solvent attraction has been recog¬ 
nized, for instance, in Ref. 5. While, in contrast to the 
ASHS model, the SAO model does not allow for an an¬ 
alytical solution in the PY approximation, its associated 
effective solute-solute pair potential can be exactly de¬ 
rived in the semi-grand-canonical ensemble, analogously 
to the case of the pure AO model.Moreover, and 
also in analogy with the AO model, ^^’^9-21 ^ pg^jj. 

potential turns out to be the only one contributing to 


3 


AHS 

ASHS 


ILO 



K 

< 

NASHS 


Z 



AO 

SAO 


0 

1 

2 


FIG. 3. Plane aaajGs vs sketching different models 
mentioned in the text. The general class of NASHS mod¬ 
els include, as limiting cases, ASHS (o-ss/ns = 1), SAG 
(osalos = 0), and NAHS ~ 0)- The intersection of the 
NAHS line with the ASHS and SAO lines define the AHS and 
AO models, respectively. In this paper, we will be concerned 
with the ASHS and SAO models. 


the exact effective interaction among the solutes if the 
size ratio q = aajai is smaller than the threshold value 
go = 2/-\/3 — 1 ~ 0.1547. A careful comparison between 
the results of the two models (ASHS and SAO) allows us 
to pave the way for an improved theoretical understand¬ 
ing of the above experiments. 

It is interesting to observe that, when the solute- 
solvent adhesion is set to zero, the model ASHS re¬ 
duces to a size-asymmetric additive HS (AHS) binary 
mixture, while the SAO model becomes the original AO 
model, these two mixtures having quite different critical 
behaviors upon varying metastable fluid- 

fluid demixing coexistence, responsible for the broad¬ 
ening at 77s >0 of the stable fluid-solid coexistence 
(0.492 < r]i < 0.543) for pure HSs {rjs —>■ 0),^^’^"^ re¬ 
mains always metastable and exists at small enough q in 
the AHS case, whereas it becomes stable at large g in the 
AO case, where a triple point appears. Figure 3 sketches 
(in the plane cTss/o's vs t^^) the different models referred 
to above. 

The organization of this paper is as follows. Section 
H presents the problem of the effective solute interaction 
mediated by the solvent particles within a general frame¬ 
work. This is followed by Sec. HI, where the PY solu¬ 
tion for the ASHS model is exploited to find the spinodal 
curves of the original mixture and the effective solute- 
solute pair potential. The exact derivation of the effec¬ 
tive potential in the SAO model with a size ratio g < go is 
addressed in Sec. IV, its three contributions being clearly 
identified. Next, the different scenarios for criticality in 
the ASHS and SAO effective systems are analyzed via the 
second virial coefficient and the Noro-Frenkel criterion^® 
in Sec. V. A more detailed analysis for the SAO model 


is performed via a first-order perturbation theory in Sec. 
VI. Finally, our findings are discussed and put in perspec¬ 
tive in Sec. VH. The most technical details are relegated 
to four appendixes. 


II. GENERAL FRAMEWORK 

Consider a colloidal binary mixture of Ns small (sol¬ 
vent) and Ni large (solute) particles, identified by the co¬ 
ordinates ..., and \ ..., r^^}, re¬ 

spectively, in a volume V. 

Assuming pair interactions, i.e., assuming the particles 
are nondeformable, nonpolarizable, ..., (see Ref. 26 for 
a recent discussion on the reliability of this assumption), 
the total potential U can be written asU = Uss+Uu+Usi, 
where 


Na-l Na 

^ -r^®^l), (2.1) 

2=1 

Ni -1 Ni 

^ 7 ^n(|rf^ - r®|), (2.2) 

2=1 

Ns Ni 

(2.3) 

i=i j=i 

The canonical free energy F{Ns, Ni,V,T) is then given 
by 


e 


-fiF 


A-3iV.A-3^< 

NslNil 


dr 


Ns 


fj^Y^ie-P(Uss+Uii+Usi)^ 


(2.4) 

where = l/ksT [ks being the Boltzmann constant), 
Ag and A; are the de Broglie thermal wavelengths asso¬ 
ciated with the small and large particles, respectively, 
and we have used the short-hand notation dr^’^ = 
with a = s, /. 

Following standard prescriptions, 3^47.27 

principle trace out all the microscopic degrees of free¬ 
dom associated with the solvent particles and recast 
Eq. (2.4) in a form of an effective one-component sys¬ 
tem for only the solute particles with a potential energy 
, r 2 \ ..., ). More specifically. 


e-/3Vf 


VglAf^ 


drN^e-P^Uss+UsO^ 


(2.5) 


so that, Eq. (2.4) becomes 

[ dr^‘e-^^‘N ( 2 . 6 ) 

In general, however, the effective potential C/))® is 
not pairwise additive, meaning that apart from pair- 
interaction terms (and less relevant zero- and one- 
body terms), it requires three-body, four-body, ... terms. 
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Thus, the general structure of Uff would be 

Ni Ni 

uf = iV,4“’ + E 4‘’ (r''>) + E 4“’(|r<'> - rf I) 

i=l i<j 

Ni 

+ ^ + (2.7) 

i<j<k 

The physically most relevant contribution is expected to 
be the one associated with the effective pair potential 

(o') 

vii{r) = v\i (r), in which case one can approximately 
neglect with n > 3. 

Now we specialize to the general class of NASHS mod¬ 
els described in Sec. I. The Pasif) and <pu(r) pair inter¬ 
actions are of HS type characterized by diameters CTss 
and ail, respectively, while the small-large interaction 
(fisiir) is of SHS type^’^® with a hard-core distance a si 
and a stickiness parameter , the latter measuring the 
strength of surface adhesiveness. Therefore, the relevant 
Mayer functions fa-y{r) = _ i are 

fssir) = -Q{ass - r), (2.8) 

fuir) =-<d{aii-r), (2.9) 

fsi{r) = -0(tTsi - r) -I- -^^5{r - Osi). (2.10) 

12Tsj 

Here, Q{x) is the Heaviside step function and S{x) is the 
Dirac delta function. To simplify the notation, we adopt 
the viewpoint of the large spheres by calling ai = an their 
diameter and defining as as the diameter of the small 
spheres as felt by the large ones, so that a si = {as+ai)/2. 
Thus, the size asymmetry of the mixture (again from the 
viewpoint of the solute particles) is measured by the ratio 
q = as!ai < 1, while the nonadditivity of the hard-core 
interactions is measured by the ratio agsjas < 1 (where, 
as said before, we discard here the case of negative non¬ 
additivity). For later use, let us introduce the partial 
packing fraction of species a as rja = TrpctCr^/b, where 
Xa = Na/N is the concentration of species a = s, I and 
Pa = Na/V is its density. The total number of particles 
and number density of the fluid mixture are N = Ni +Ns 
and p = N/V, respectively. 

As discussed in Sec. I, we now particularize to two in¬ 
teresting particular cases that are identified by the ratio 
ass/as', the ASHS model (where ass/as = 1) and the 
SAO model (where o'ss/o's = 0). The first model was 
studied before by two of us (it was called System A in 
Sec. V of Ref. 12) and has been rejuvenated by a recent 
study by Chen et al.^^ The second model is an exten¬ 
sion of the well-known AO model, except that a sticky 
(or adhesive) interaction exists between the solvent and 
the solute particles. To the best of our knowledge, it has 
not been studied before. In both cases we will be able to 
derive the effective pair potential vu (r) = vp ' (r) [see Eq. 
(2.7)] either within the PY approximation in the canon¬ 
ical ensemble (ASHS model) or in an exact way in the 
semi-grand-canonical ensemble (SAO model). 


III. THE PY APPROXIMATE SOLUTION OF THE 
ASHS MODEL 

The solution of the PY approximation for the ASHS 
model was recently studied in Ref. 12. The PY solution 
actually extends to the more general formulation where 
the Baxter stickiness coefficient^’^® between a particle of 
species a and one of species 7 is Since here 

we choose Tss —>■ 00 and tu —> 00 , we can only have 
adhesion between unlike particles and >0 measures 
its strength. 


A. Spinodal curve 


From Eq. (85) of Ref. 12 we find the following expres¬ 
sion for the spinodal of the full binary mixture in the 
(jJstVi) plane, as obtained from the PY approximation. 




1 + (1 + g)(i -Vs- vi)/3{ps + qyi) _ ^ 

- 1 +^( 1 +^) ( 1 +’^) - 


{^ + q){i 1 s + qyi) 
mi - Ps - vi) 


(3.1) 


which, as it should, is symmetric under the exchange 
Vs ^ Vi and q O 1/q. For a fixed q, there is a maxi¬ 
mum value of m/ for which Eq. (3.1) admits a solution 
with ? 7 s > 0 and pi > 0. We will denote this maximum 
value with fsi and the corresponding solution, the criti¬ 
cal point, with (fjs,vi). In particular, at g = go we find 
Tsi = 0.014448, ps = 0.019 839, and pi = 0.101645. For 
Tsi < Tsi the solution of Eq. (3.1) is a closed curve in 
the {ps, pi) plane within which the thermodynamically 
unstable region lies, as shown in Fig. 4. As we can see, 
the spinodal curve does not change much for Tsi < 0 . 001 , 
where it is crossed by the straight line representing the 
critical packing fraction ps = Ps- These findings are in 
complete agreement with those reported in Ref. 11. 

Note that Eq. (3.1) is a particular case of an equation 
for a general mixture derived by Barboy and Tenne,^® 
that should however be handled with great care.®° 


B. Approximate effective one-component fluid 

As explained in Sec. H, one could in principle inte¬ 
grate out the solvent degrees of freedom to obtain the 
effective solute potential C/f® [see Eqs. (2.5) and (2.7)]. 
Here we want to focus on the pair interaction poten- 
tial vii{r) = v/ (r). This function can be identified 
from the solute-solute radial distribution function gii{r) 
in the infinite dilution limit {xi —>■ 0 ) since in that limit 
only pair interactions contribute to gii{r). Therefore, 
guir) —>■ and hence 

mm) = - lim \ngiiir). 

xi->0 


(3.2) 
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FIG. 4. PY spinodal for q = qo and several values of Tgi. The 
straight line is rja = qt and the circle is the critical point at 
Tsi = 0.014 448. 


In the limit of no adhesion [tsi —^ oo), viiir) becomes the 
usual depletion potential.For further use, we will 
refer to entropic regime as the one with Tgi ^ 1, close to 
a size-asymmetric binary HS mixture. Reciprocally, the 
nonentropic regime will refer to a system with a small 
Tg/. The transitional regime will correspond to Tgj ~ 1. 

Since 7]s is supposed to be finite in Eq. (3.2), it is not 
possible to obtain the exact effective pair potential viilr). 
On the other hand, it can be obtained again from the 
PY solution, as described in Appendix B. Note that, al¬ 
though the infinite dilution limit is applied as a short¬ 
cut to derive the pair potential u/;(r), at a nonzero so¬ 
lute concentration the full effective many-body poten¬ 
tial Ufj® includes nonpairwise terms, as represented by 
v\i ’{Yi,Vj,Vk) and higher-order terms in Eq. (2.7). 

In Eig. 5 we report a few representative examples of 
the effective solute-solute pair potential corresponding 
to the ASHS model in the PY approximation (see Ap¬ 
pendix B). Figure 5(a) shows the influence of the solute- 
solvent stickiness at fixed ps = 0.1 and q = qq. One 
can clearly observe the different shape of the potential 
in the entropic {tsi = 10^), transitional {tsi = 1), and 
nonentropic {tsi = 0.12) regimes. In the former case 
{tsi = 10^), the potential is essentially attractive (ex¬ 
cept for a slight hump in the region rjoi < 1 -|- g), thus 
reflecting the depletion mechanisms (see Fig. 1). More¬ 
over, at this very high value of Tg;, the discontinuity of 
the potential at rjai = 1 + q [see Eq. (B15)] is not vis¬ 
ible. In the transitional regime (rgj = 1), however, the 
discontinuity at r/cr; = 1 -I- g is already noticeable and 
the potential in most of the inner region 1 < rjai < 1 -fg 
has changed from attractive to repulsive. These two 
features are widely enhanced in the nonentropic regime 
{tsI = 0.12): there is a high discontinuity at r/ai = l + q 
and the effective potential is strongly repulsive in the 
whole region 1 < rjoi < 1 -I- g. Eurthermore, a strong 
repulsion appears as well in the outer region rjoi > l-fg. 
Figure 5(b) shows that an increase of the solvent density 
magnifies the characteristic features of the effective po- 



FIG. 5. Effective solute-solute pair potential in the ASHS 
model, as obtained from the PY approximation (see Appendix 
B). In panel (a) the stickiness parameter is varied at a fixed 
solvent packing fraction pa = 0.1, while in panel (b), pa is 
varied at fixed Tai = 0.12. In all the cases the size ratio is 
q = go. 


tential in the nonentropic regime. The physical origin 
of the repulsive regions in the nonentropic regime can 
be ascribed to the net pulling role played by the solvent 
particles attached to the two solutes. This effect will be 
identified more clearly in the SAO model (see Sec. IV). 
As for the (attractive) discontinuity at rjai = 1 -|- g, it 
can be attributed to the bridging effect of solvent parti¬ 
cles attached to both solutes. This bridging mechanism 
is absent if r/ai = (1-l-g)''" but appears if r/ai = (l-|-g)“. 

Dijkstra et al.^®’^^ already showed that the effective 
potential in the entropic regime is unable to produce a 
stable demixing phase transition with reasonably small 
g. On the other hand, the step attraction at r/cr; = 1 -|- g 
in the potential associated with the nonentropic regime 
can lead to a demixing transition, as shown in Ref. 12. 
This is the phase instability studied in the {ps,pi) plane 
in Sec. Ill A. 

With all due cares, the shape of the effective potential 
in the nonentropic regime depicted in Fig. 5(b) can be 
schematically represented as a square-well (SW) poten¬ 
tial of width qai and depth e ^ |u;/(cr;(l -I- g)“)|, with 
an additional repulsive tail starting at r = cr;(l -|- g)^. 
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We can then exploit the fact that the phase behavior 
of a one-component SW fluid is well established. 

For example, it is sufficient to heuristically consider 
the approximate critical value^® of the reduced tem¬ 
perature T* = ksT/e to find the appearance of an 
open phase coexistence region at high rj^ (well separated 
from the closed one predicted in Sec. Ill A at low rjs). 
This coexistence region is known to be present in the 
highly asymmetric AHS mixture^®’^^ (be., for small q 
in the limit —>■ oo). The effective problem proce¬ 

dure that we followed suggests that, quite intuitively, 
such a region will not disappear when the attraction 
is switched on at small Tgi. It is interesting to ob¬ 
serve that such a reentrance at large rjs is not predicted 
by an analysis of the behavior of the effective second 
virial coefficient [see Eq. (B16)], according to which 
l/T* = \n[l + {l-Bf/Bf^)/{3q + 3q^+ q^)], where 
SHS = 

^crf. The two heuristic criteria based on an 
effective SW temperature T* agree quite well for small 
values of rjs [as expected from the curve rjs = 0.01 in Fig. 
5(b)], but the Bf^ criterion presents a diverging T* at 
a value of rjs such that = B^^ and becomes mean¬ 
ingless thereafter (i.e., when > B^^). For instance, 
ii q = qo and Tsi = 0 . 12 , the condition > B^^ is 
satisfied for rjg > 0.274. The fact that B^^ > B^^ if r/g 
is large enough is directly related to the increase of the 
effective size of the dressed solute particles, as depicted 
in Figs. 2(c) and 2(d). 

Of course, the effective one-component fluid is not fully 
equivalent to the original binary mixture, as we are ne¬ 
glecting three-body (and higher) terms in the effective 
total potential [see Eq. (2.7)]. Moreover, the potentials 
of Eig. 5 are the outcome of the PY approximation. Yet, 
they are expected to give reasonable approximate results 
in the spirit of an effective fluid. Chen et al.^^ devised a 
similar approximate mapping of the PY solution for the 
true binary mixture onto a one-component SHS model, 
from which they were able to read-off the binodal us¬ 
ing accurate Monte Carlo (MC) results by Miller and 
Erenkel.^^ 

While some caution must be exercised when using the 
pairwise potential formally obtained in the limit 77 ; —>■ 0 
to predict the phase diagram at finite rji , this keeps being 
a useful procedure to reduce the complexity of the binary 
mixture problem,^® allowing one to get additional physi¬ 
cal insight without the need, for example, of performing 
computer simulations of the full binary mixture. 


IV. THE SAO MODEL 

As shown in Sec. Ill, the ASHS model (ctss/cs = 1) 
admits a PY analytical solution but only an approxi¬ 
mate reduction to an effective one-component fluid. The 
SAO model (uss/o's = 0) is, in some sense, complemen¬ 
tary to it, as it does not admit an analytical solution, 
not even in the PY approximation, but it does admit an 
exact reduction to an effective one-component fluid for 


q < qo = 2/-\/3 — 1 — 0.1547, when a solvent particle 
can fit into the inner volume created by three solutes 
at contact,®® so that a solvent particle cannot overlap 
simultaneously with more than two (nonoverlapping) so¬ 
lute particles. This corresponds to q < 1 in one spatial 
dimension.^® 

To proceed, it is convenient to change from the canon¬ 
ical {Ns, Ni,V,T) ensemble to the semi-grand-canonical 
(/is, Ni, V, T) ensemble,where jig is the chemical po¬ 
tential of the solvent component. The corresponding 
thermodynamic potential B{jis,Ni,V,T) is constructed 
via the Legendre transform 

B{jis,Ni,V,T) = F {{Ns),Ni,V,T)- jig{Ng). (4.1) 
Thus, the counterpart of canonical Eq. (2.4) is 


-/3.F 


= E 


.Ns 


Ns=0 

1 

TVilAp' 


Ng\Ni\Af^' 




where 




■Xf 


(4.2) 


(4.3) 


is the solvent fugacity and 


^Ns 


dr 

Ns=0 


Ns^-PUsi_ 


(4.4) 


Note that in Eq. (4.2) we have taken into account that 
Ugg = 0 in the SAO model. 

Inserting Eq. (2.3) into Eq. (4.4) it is easy to obtain®^ 



where 0 = — Uu represents the grand potential of 

an ideal gas of solvent particles in the external field of a 
fixed configuration of Ni solute particles with coordinates 
{rf^}. Expanding in products of Mayer functions, can 
be written as 


^max (<?) 

n= (4-6) 

n—0 

Here, is the contribution to H stemming from the 
product of n Mayer functions fgi. The upper limit 
’T-max('?) is the maximum number of nonoverlapping so¬ 
lutes that can simultaneously overlap with a single sol¬ 
vent particle. For n > nmax(9), at least one of the factors 
fsi vanishes and so does n„. If 9 < qo, then nniax(<z) = 2, 
implying that the exact effective potential C/fj® does not 
include three-body (or higher order) terms. In the inter¬ 
val <70 < 9 < Ij ?^max(9) grows by steps as q increases, 
reaching a maximum value Umax ( 9 ) = 11 (since a solvent 
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r < <7i, 

ai <r < (Tj(l + q), (4.12) 

r > ai{l + q), 



FIG. 6. Cartoon describing the three effects (depletion, 
pulling, and bridging) contributing to the effective solute- 
solute interaction in the SAG model. The grey large spheres 
represent the solutes of diameter ai at a distance r < ai + 
Us= -I-q). They are surrounded by a sea of smaller 
spheres (the solvent) of diameter Us = gcr; {q = 0.2 in the car¬ 
toon) that can overlap among themselves and have a sticky 
surface (represented by a thick perimeter). Some of the sol¬ 
vent particles (the yellow ones) do not touch the solutes and 
so they contribute to the (attractive) depletion effect, which 
is represented by V'd(r), a volumetric term independent of 
Tai. Other solvent particles (the orange ones) are adhered to 
one of the big spheres, thus contributing to the (repulsive) 
pulling effect, represented by tpp{r), which is a surface term 
proportional to . Finally, some other small particles (the 
red ones) are adhered to both solutes, giving rise to the (at¬ 
tractive) bridging effect, represented by 'iph{r), which is a line 
(intersection of two surfaces) term proportional to 


particle can simultaneously overlap with 12 nonoverlap¬ 
ping solutes only if q > 1). The hrst few terms in Eq. 
(4.6) are 


/30o =-ZsV, 


/30i = -Zs 



drfsi{\r - 




I), 


/ 3 O 2 = -2s 





\)fsl{\ 


r — 



(4.7) 

(4.8) 


N, 



<^„(|rl'^-rW|) (4.9) 


Equation (4.9) allows us to identify the exact effective 
pair potential as 


l3vii{r) = P'fii{r) - Zs J drs fsi{rs)fsi{\rs - r\). (4.10) 
Now, making use of Eq. (2.10), one can obtain 


{ 00, 
V’(r-), 
0 , 


where = Zs{TT/6)ag is the (nominal) solvent pack¬ 
ing fraction of a reservoir made of noninteracting solvent 
particles and 


ij){r) = tpdir) + tjjp{r) + ijjb{r) 


(4.13) 


with 


fpd{r) = -^ [ dvs <d{asi - rs)Q{asi - |rs - r|) 

TTCrf J 

{l + q- rla{f{2 + 2q + r/ai) 


2g3 


(4.14) 


V'p(»') = [ dVs S{rs - asi)Oicrsi - |rs - r|) 

{l + qf{l+q-r/ai) 


dg^Ts, 


(4.15) 


2 p 

i’hir) = - 24^^3.^2 J “ crsz)( 5 (|rs - r| - Osz) 


a + dY 


192q^T^ir/ai' 


(4.16) 


The effective solute-solute force fu{r) = —dvu{r)/dr 
(outside the hard core, r > oi) is 

= - [V'd(0 + V'p(0 + 0(1 + g - r/ai) 


Vs 


(1 q)^ 


( 4 . 17 ) 


where the delta term reflects the discontinuity of vii{r) 
at r = tTi(l -|- q) and 


Yd(.r) = 


2q^ai 


(1 + gf- 


'I J 


(4.18) 


= . (4.19) 


4q^aiTsi ’ 


192g3T-2 r^/cr;' 


If <7 < goj the general relationship between the reservoir 
packing fraction vY'^ (or, equivalently, the fugacity Zs) 
and the values Vs and vi of the binary mixture is derived 
in Appendix C with the result 


(r) 

Vs = Vs 


i-.„(i+<,)qi-jE)-l^ 


t■(Tl{l+q) 

X / drr‘^il){r)gef[{r\vuvYY 

J ai 


(4.20) 


jdVti 


ZsViV{l + qY 



(4.11) 


where gesiflVhvYY is the radial distribution function of 
a pure fluid of large particles interacting via the effective 




















r/0[ 



r/Oj 



r/0[ 



r/Oj 


FIG. 7. Plot of 'ip{r) = I3vu(r)/T]i''* for (a) t^i = 1, (b) Tsi = 0.2, (c) Tai = 0.025, and (d) Tsi — 0.02. The solid lines correspond 
to the SAG model at the threshold value q = qo, while the dashed lines correspond to the SWAO model with q = qo — 
Asi = ^/{l + q), ^ = 10~^. The insets in panels (c) and (d) show magnified views of the curves for r/ai < 1 + go- 



FIG. 8. Plane Tsi vs r showing the region where the effective 
pair potential in the SAG model for the threshold value g = go 
takes positive values. Gutside the shaded region the potential 
is negative. 


pair potential vu (r) at a packing fraction rji. Up to second 


order in rji, Eq. (4.20) becomes 


Vs ~ V> 


(r) 


1 - vii^ + qf 1 - 


4 ts, 


,2^3 


I2r;/g' 


f-rjiil+q) 


dr r'^'ip{r)e 


(4.21) 


Interestingly, exact effective pair-potential (4.12) can 
be equivalently obtained from a density expansion of the 
approximate PY effective potential of the ASHS model 
described in Sec. Ill, upon neglecting terms of order 

Cr') 

higher than linear in rjs and identifying rjs with 77 ) , 
that is correct in the solute infinite dilution limit [see 
Eq. (4.20)]. This is not a coincidence^^ because the PY 
theory gives the exact radial distribution function to first 
order in density (and therefore it gives the exact effective 
potential to that order) and the relevant Mayer diagram, 
containing only one solvent particle, is the same whether 
the mixture is additive or not. 

The three terms appearing in Eq. (4.13) bear a par¬ 
ticularly simple and instructive physical interpretation. 
The first term, tpdir), [see Eqs. (4.14) and (4.18)] is the 
conventional AO effective potential.^® If r < cr;(l -|- g), 
no solvent particles fit in the line joining the centers of 
the two solute particles. This is the typical configuration 
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of depletion when the solute-solvent interactions are of 
HS type, giving rise to an effective attraction between 
the solutes (with a force decreasing its strength quadrat- 
ically with increasing distance). Now imagine we switch 
the stickiness on. Interestingly, this produces two com¬ 
peting effects. Firstly, the solvent particles attached to 
the outer surfaces of each facing solute tend to pull the 
solutes apart, producing an effective solute-solute repul¬ 
sion with a constant force strength. This is represented 
by the “pulling” term [see Eqs. (4.15) and (4.19)]. 
Secondly, the solvent particles attached to both facing 
solutes (the “bridges”) tend to increase attraction (with 
a Coulomb-like force strength decreasing with increasing 
distance), this bridging effect being represented by the 
term 'il}h{r) [see Eqs. (4.16) and (4.19)]. These three ef¬ 
fects are schematically synthesized in Fig. 6. 

It is interesting to remark that the SAO model can be 
easily extended by replacing the solute-solvent sticky sur¬ 
face by a finite-width (Ag;) SW interaction. The result¬ 
ing SWAO model is worked out in Appendix D. In this 
case, the condition for an exact reduction of the effective 
solute interaction to pairwise terms is 9 ( 1 - 1 - Agj) -|- Ag/ < 
90- 

The interplay of the three contributions to ^/’(r) gives 
rise to interesting transitions in the shape of the depletion 
potential, as illustrated in Fig. 7 for the SAO and SWAO 
models. Let us comment the curves corresponding to the 
SAO model. For relatively weak stickiness, < tZ^ = 
24(1 — yjl — q — I2)I{1 9 )^, the pulling effect domi¬ 

nates over the bridging effect for all distances but is dom¬ 
inated by the depletion effect, except for distances close 
to r = ai{lq). Consequently, the effective potential is 
attractive near r = ui and repulsive near r = cri(l-|- 9 ), as 
happens in Fig. 7(a). Next, in the intermediate regime 
tZ^ < < T+^ = 24(1 -H a/1 - 9 - 92 / 2 ) 7(1 -f qf the 

pulling effect dominates for all distances and the poten¬ 
tial is purely repulsive, except for the discontinuous jump 
at r = cr/l + 9). This is represented by the case of Fig. 
7(b). In the strong stickiness regime < 48 the 

depletion effect is practically irrelevant and the pulling 
effect is dominated by the bridging one, except in the 
region r < ai{l q). As a consequence, the effective 
potential is slightly attractive near r = ai and slightly 
repulsive near r = ct /1 -I- 9), as happens in Fig. 7(c). Fi¬ 
nally, for very strong stickiness (r/)^ > 48) the bridging 
dominates over the pulling for all distances and the po¬ 
tential is purely attractive. This is the case displayed 
in Fig. 7(d). Those features are essentially preserved 
in the case of the SWAO model, except that the jump 
at r = aiil -\- q) is replaced by a rapid (but continu¬ 
ous) increase of the potential between r = ai{l q) and 
r = cr/l -I- 9)(1 -I- Agj). 

From Eqs. (4.14)-(4.19) it is easy to see that in the 
SAO model the effective potential and force are positive 
if 


6r — W6r(5f — 2) 6r -|- W6r(5f — 2) , ^ 

2Af{l-f){2 + f) - - 24f(l-r)(2-hf) ^ ^ 


and 


2f - v^2(3r2 - 1) 2f + v^2(3r2 - 1) 

24r(l — 72 ) “ “ 24r(l — f2) ’ 


(4.23) 


respectively, where f = r/[CT/l -I- 9 )]. Figure 8 shows 
the region in the plane Tsi vs r where ^/(r) >0 for the 
threshold value 9 = 90 . 

As can be seen from Figs. 7(c) and 7(d), the effec¬ 
tive potential in the regime of strong stickiness clearly 
resembles that of a SW potential of width qai and depth 
/3e = 77 /^|V'(o'i(l -I- 9 ))| = 77//1 -I- 9“/3/192t/,. In this 
case, the effective phase behavior of an equivalent SW 
fluid would suggest that, for a given 9 < 90 and a suffi¬ 
ciently small Tsi , we have the appearance of just one lower 
critical point ( 77 /^'^, 77 /) and the instability region does 

not close itself again at 77 /^ > 77 /^'^. This is a scenario 
quite different from the one in the model ASHS, where 
we found at least one closed island with a lower and an 
upper critical point (see Sec. Ill A). It would then be suf¬ 
ficient to switch on a hard-core repulsion (with ass = 0 's) 
among the solvent particles to have a closed spinodal. 
Along similar lines, it is also interesting to observe that 
the threshold packing fraction rj* defined in Sec. I clearly 
diverges in the SAO model because the solvent particles 
can freely overlap. 


V. THE NORO-FRENKEL CRITICALITY CRITERION 


In 2000, Noro and Frenkel (NF)^® argued that the re¬ 
duced second virial coefficient B 2 /B^^, rather than the 
range and the strength of the attractive interactions, 
could be the most convenient quantity to estimate the 
location of the critical point for a wealth of different col¬ 
loidal suspensions. Their criticality criterion for particles 
with variable range attractions,^® complemented by the 
simulation value of the critical temperature obtained in 
Ref. 37 for the SHS model, yields B 2 /Bf^ ~ —1.21. 

In this section we apply the NF criterion to the two 
models discussed before: the ASHS model (see Sec. Ill) 
and the SAO model (see Sec. IV). In both cases, if vii{r) 
is the effective solute-solute pair potential, the associated 
second virial coefficient B^^ is given by 


Bf 

RHS 


pOO 

/ dr 

^-0vii{r) _ ^ 




(5.1) 


where B^^ = 27rcr//3 is the virial coefficient for HSs of 
diameter < 7 /. Paradoxically, while the explicit PY expres¬ 
sion of j3vii {r) in the ASHS model is rather cumbersome 
(see Appendix B), its associated second virial coefficient 
is much easier to obtain thanks to properties of the 
Laplace representation. The result can be found in Eq. 
(B16). In contrast, in the SAO model the exact expres¬ 
sion of Pvii{r) is very simple [see Eqs. (4.12)-(4.16)] but 
the computation of B^^ needs to be done numerically. 
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FIG. 9. Second effective virial coefficient as a function of for q — qo and several values of for (a) ASMS and (b) SAO 
models. The thick horizontal line corresponds to the NF criticality criterion = —1.21. Panels (c) and (d) show the 

dependence of the critical value 77 ^’'^ = (according to the NF criterion) as a function of Tsi for the ASHS and SAO models, 
respectively. Note that a logarithmic scale is used on the abscissas. 


It is particularly instructive to observe that the NF 
criterion confirms the very different critical behavior be¬ 
tween the ASHS model {(Jss/o-s = 1) and the SAO model 
{ass/ 0 's = 0). In Figs. 9(a) and 9(b) we compare the 
second effective virial coefficients for the two models as 
functions of Tsi for q = qo and several values of 77 s ■ Here 
we have identified rjs in the ASHS case, in consis¬ 

tency with the fact that the effective potential is derived 
in the infinite solute dilution limit. The loci of points in 
the plane vs Tsi where = —1.21 are dis¬ 

played in Figs. 9(c) and 9(d). Inside the shaded regions 
one has < —1.21 and thus phase coexistence is 

possible, according to the NF criterion. 

As we already knew from the residts of Sec. HI A, Figs. 
9(a) and 9(c) show criticality in the ASHS model only for 

sufficiently small Tsi and qs ■ On the other hand, the sce¬ 
nario present in the SAO model is completely different. It 
is easy to check that a critical point in the pure AO model 
{tsi —>■ 00 ) exists only, according to the NF criterion, if 
> 0.318. However, the presence of stickiness (finite 
Tsi) dramatically changes the picture. For any q) , there 
exists a critical point if Tsi is small enough. Beyond a cer¬ 
tain threshold value, criticality abruptly disappears and 
then (only if qi^'^ > 0.318) it re-enters at a sufficiently 


large value of Tsi- Thus, if q^s'^ > 0.318 there exists a win¬ 
dow of values of TsI where no phase separation is possible. 
Note that values of q^s'’ > 1, as displayed in Fig. 9(d), 
are not unphysical in the SAO model since the reservoir 
consists in an ideal gas of noninteracting small particles. 

It must be remarked that the bridging and pulling ef¬ 
fects are more important in the nonadditive SAO case 
than in the additive ASHS one, since in the latter the 
mutual exclusion of solvent particles interferes with their 
ability to attach to the solutes. As illustrated in Fig. 9, 
this leads to paramount differences in the critical behav¬ 
ior of the two extreme models. For intermediate NASHS 
models with 0 < ass/og < 1 (see Fig. 3) a transition 
from Figs. 9(b) and 9(d) to Figs. 9(a) and 9(c), respec¬ 
tively, can be expected as the excluded volume among 
the solvent spheres is gradually increased. 

Note also that in the ASHS model [Fig. 9(a)] the results 
are approximate (PY) and the solute concentration is 
zero. On the contrary, in the SAO model [Fig. 9(b)], 
the results are exact and valid for any finite solute and 
solvent concentrations. While both models coincide in 
the limit of vanishing solvent concentration, in practice 
this equivalence requires extremely small values of 77 ) . 
For instance, at q^s'^ = 10“® both values of differ by 
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nearly 2%. 


VI. PERTURBATION THEORY FOR THE SAO MODEL 


From Sec. V we conclude that the “hidden” fluid-fluid 
phase separation observed by Dijkstra et al.^® in their 
study of the AO model could be stabilized by adding ad¬ 
hesion, as in our SAO model. This can be quantified more 
precisely using a first-order thermodynamic perturbation 
theory.^® 

Assuming the HS fluid as reference system, we can 
write the Helmholtz free energy per particle of the effec¬ 
tive solute system as 

/3Fhs , ,, M p(i+9) 

-^ = —^ + drr 

( 6 . 1 ) 

where PF^s/Ni = {4r]i - Srjf)/{1 - rji^ + \n{r]i) + const is 
the Carnahan-Starling^^ HS expression, is given by 
Eqs. (4.13)-(4.16), and ^hs is the HS radial distribution 
function in the PY approximation"*^^, which in the interval 
ai < r < ai{l + q) < 2 can be written as 


9m{r\r]i) 


y lim 


[t - tn{r]i)]tL{t\r]i) ed’’ P 
S{t\r]i) r 


( 6 . 2 ) 


where we are measuring lengths in units of ai, 

S(t\r]i) = (1 - -f 6r]i{l - r]i)f + ISgft 

-I2iqiil + 2r,i), (6.3) 

L{t\rii) = {l + rii/2)t + l + 2r]i, (6.4) 


and tn{r]i) {n = 1,2,3) are the zeros of S{t\rii). The 
hrst-order Helmholtz free energy of Eq. (6.1) can thus be 
calculated analytically. 

The compressibility factor Z = f3p/p and chemical po¬ 
tential /i are then found through 


d{PF^^/Ni) 

^ = TJl 


= Z®” -h 


■eff 


dm 


(r) 

Vs 


(6.5) 

( 6 . 6 ) 


The critical point pf) is determined by numerically 

solving the following set of equations: 


d{mz<^^) 


dm 


= 0 , 


(^)c c 

Vs ,Vi 


d^mz^^) 

dvf 



= 0 . 


(6.7) 

( 6 . 8 ) 


In Fig. 10 we show the critical point , vf) for the 
fluid-fluid coexistence in the SAO model at the threshold 
value q = Qq as SL function of Tgi ■ The figure confirms the 
scenario predicted in Sec. V from the NF criterion. In 
fact. Fig. 10 shows a relevant mutual consistency between 



FIG. 10. Critical point for the fluid-fluid coexistence in the 
SAO model for g = go as a function of TsI ■ The lines with sym¬ 
bols are obtained from perturbation theory, while the solid 
line corresponds to the NF criterion [see Fig. 9(d)]. A loga¬ 
rithmic scale is used on the abscissa. Equation (4.21) is used 
for the conversion between the reservoir and the solvent den¬ 
sities. 


the curves for as obtained from both independent 
approaches. There is a range of adhesion for which there 
is no criticality. For high adhesions (small Tsi) we have 

fr) 

phase coexistence in the region of low rjs region, while 
for low adhesions (large t^i) the criticality exists in the 

(r) 

region of high rjs ■ Of course, we expect a breakdown of 
the perturbation theory treatment as soon as stickiness 
becomes too strong. Also, as soon as g > go we are 
neglecting three-body (and higher) terms. 


VII. CONCLUSIONS 


In this paper, we have studied two complementary 
models of a binary mixture of (small) solvent and (large) 
solute particles, where in both cases unlike particles ex¬ 
perience an attractive adhesion interaction of Baxter’s 
type.^ We studied the derivation of an effective solute- 
solute pair-potential for the two models in the regime 
of large size asymmetry (g = CTs/ct; <C 1) and discussed 
analogies and differences of the corresponding phase be¬ 
haviors, as obtained from the resulting effective one- 
component fluid. 

In the first model, that we dubbed ASHS, both solute- 
solute and solvent-solvent particles interact as HSs and 
the reduction to an effective one-component fluid can be 
carried out only approximately via a small solute den¬ 
sity expansion. By contrast, this model admits an exact 
analytical solution within the PY approximation. In the 
limit of vanishing solute-solvent adhesive attraction, this 
model reduces to the usual AHS binary mixture, that is 
known not to display any phase separation within the 
PY approximation. This might, however, be ascribed to 
the limitations of the PY closure, as other more sophisti¬ 
cated theories, as well as numerical simulations, support 
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the existence of phase separation, albeit metastable with 
respect to freezing, at sufficiently large concentrations 
and size asymmetry (in this context, nevertheless, see 
Ref. 43). In this case, our analysis of the ASHS model 
confirms previous findings of a similar study by Chen et 
al.^^ in predicting a closed region in the (?7s,?7i) plane 
where phase separation occurs. 

While the ASHS model has been around for some 
time,^^ the second model (denoted as SAG) is, to the 
best of our knowledge, new. In this case, solvent parti¬ 
cles behave as an ideal gas within each other —but still 
they experience a SHS interaction with the solutes. In 
the limit of no adhesion between solute and solvent, this 
model reduces to the well-known AO one, and we have 
extended the analysis performed by Dijkstra et al.^^ to 
the present case. As in the AO case, even in the SAO 
case the solvent degrees of freedom can be traced out ex¬ 
actly above a well defined size asymmetry (that is, below 
a critical value qo of the size ratio q), so that the result¬ 
ing effective one-component pair potential is exact. By 
contrast, it is not possible in this case to obtain an ex¬ 
act analytical solution of the binary problem (not even 
within the PY approximation), so we resorted to study 
a first-order thermodynamic perturbation theory of the 
corresponding exact effective solute-solute pair potential. 

In both models, effective potentials can be explained 
in terms of “pulling” and “bridging” effects in addition 
to the usual “depletion” mechanism. In the SAO case, 
the analytical expressions of the effective potential de¬ 
rived in Eqs. (4.12)-(4.16) allow for an interesting direct 
physical interpretation. The pulling effect is represented 
by the term proportional to [see Eq. (4.15)], as the 
same (solvent) particle must be in contact with one of 
the solutes and outside the exclusion volume of the other 
solute. On the other hand, the bridging effect is repre¬ 
sented by the term proportional to [see Eq. (4.16)], as 
the same (solvent) particle must be in contact with both 
solute particles. These effects are present in both mod¬ 
els, but they are more important in the SAO case than 
in the ASHS case, since in the latter, the mutual exclu¬ 
sion of solvent particles interferes with their ability to be 
attached to the solutes. In fact, the situation sketched in 
Eig. 2(d) is inhibited in the SAO model, as represented 
by Fig. 9(d), which shows always phase coexistence at 
increasing for any fixed Tsi. 

The derivation of the exact SAO effective potential 
has allowed us to clearly assess the dramatic influence 
of solute-solvent attraction on the conventional AO de¬ 
pletion potential. This complements a recent study,^^ 
where softness in the solute-solvent repulsion was seen to 
strongly enhance the depletion mechanism. 

Leaving aside the issue of the metastability with re¬ 
spect to the fluid-solid transition, the resulting pic¬ 
ture confirms the significant impact of nonadditivity on 
the fluid-fluid phase diagram, as synthesized by Fig. 9. 
Within the NF criticality criterion, the SAO model is ex¬ 
pected to display a reentrant phase transition in terms of 
Tsi, whereas the ASHS model is not. On the other hand. 


the results for the ASHS model are compatible with a 
reentrant phase transition in terms of r]s not observed in 
the SAO model. A first-order perturbation theory on the 
SAO model confirms this picture. 

Our findings nicely confirm and complement those by 
Chen et al.,^^ but extend them to encompass a direct 
connection with the AO original model, that was miss¬ 
ing in the above study, thus paving the way to a more 
direct interpretation of the experimental results reported 
in Refs. 8 and 9. 

While direct numerical simulations of binary mixtures 
with large size asymmetries are notoriously difficult, it 
would be interesting to study with numerical experi¬ 
ments whether adhesion gives rise to the appearance of 
a metastable fluid-fluid coexistence at large solvent den¬ 
sities for the ASHS model with large q and for the SAO 
model with very small q. In addition, they open a num¬ 
ber of interesting perspectives for future studies. Even 
without resorting to a direct numerical simulation calcu¬ 
lations, a number of different theoretical approaches can 
be exploited to make further progresses. 

As the attraction between the unlike spheres vanishes 
(tsj —>• oo), the PY solution of the ASHS model reduces to 
the well-known PY solution for a binary AHS mixture,"^® 
which does not show phase separation for any size ratio, 
in spite of the possible depletion interactions. As said 
above, this seems to be an artifact of the PY approxima¬ 
tion, as shown by numerical simulations of the (approx¬ 
imate) effective one-component fluid^® and by numeri¬ 
cal solutions of the Rogers and Young (RY) closure.^® 
Thus, one possibility would be to use the RY closure on 
a binary mixture with HS interactions between like par¬ 
ticles and a short-range SW attraction (in the regimes 
where this can be considered sticky- like"*’’’”"'^®) between un¬ 
like spheres. Work along these lines is in progress and will 
be reported elsewhere. 

Another possibility would be to consider a binary 
ASHS mixture with HS interactions between small 
spheres, weak SHS interactions between the large 
spheres, and stronger SHS interactions between small and 
big spheres. This two-component model (which is known 
to be free from the thermodynamic inconsistency affect¬ 
ing the one-component model^®) could be solved rather 
easily within the PY approximation, as done for instance 
by Zaccarelli et al.®® 

Finally, it would be nice to extend the study reported 
here for the ASHS and SAO models to a more general 
NASHS model where one could tune the solvent-solvent 
diameter from zero (SAO model) to the additive value 
(ASHS model), thus encompassing both models into an 
unified framework. MC simulations for a binary ASHS 
mixture have been performed by Jamnik,®^ but not for 
the determination of the phase diagram, which has been 
studied for the one-component SHS fluid by Miller and 
Frenkel.^’’’®^ To the best of our knowledge, no numerical 
experiment has ever been tried on the NASHS binary 
mixture. 
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Appendix A: A simple geometrical argument related to Fig. 
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We hrst estimate how many small spheres of diameter 
(Ts are necessary to cover the surface of a large sphere of 
diameter ai. Assuming g = Us /ct; <C 1, the small spheres 
will be distributed on the large sphere surface approx¬ 
imately with a hexagonal packing corresponding to an 
area fraction 77hex = fa 0.907. Thus, 77hex = 4>ajA, 

where (j) is the number of the adsorbed small spheres, 
a = {'K/A)a1 is the area of the projected disk associated 
with each solvent sphere, and A = is the surface area 
of the solute particle. Therefore, 

A 27r _2 , , 

The critical volume fraction 77 * at which all large col¬ 
loidal spheres, distributed with a volume fraction rji = 
{Tr/6)piaf, can be covered is 

v*s = (A2) 


LfJ = Lf) = - -+ 


377s 


1 - 77s (1 - 77s)2 ’ 


(B7) 


^ 1 - (677s/g^)T^f 377 s 

“ l-77s 2(7(1-?7 s)2’ 


(B 8 ) 


.(1) ^ 1 + g - (1277s/g)L^f 377s 

2(1 -77s) ^2(1-77 s)2’ 


(B9) 


( 1 ) _ 1 + g 


LY> = 


377s 


2(1 -77s) 2(1-77 s)2' 


(BIO) 


(1) ^ 1 + g/2 

^( 1 - 77 ) 2 ’ 


(Bll) 


and this leads to the expression reported in Sec. I. 


Appendix B: Analytical PY expressions for the ASMS model 
in the limit xi ^ 0 

The Rational-Function Approximation (REA) 
methodology^®’^^’®^ is known to give access to an¬ 
alytical formulae of the PY solution for the ASHS 
model. 28 .j-jjjg Appendix we assume the infinite 
dilution limit for the solutes {xi —>■ 0 ). 

According to Eq. (36) of Ref. 31, the Laplace transform 
Giiis) = dr e~'‘^rgu{r) of rgii{r) is, in the limit xi 
0 , 


Gu{s) = 


Lii{s) + Lis(s)- 


Asi{s) 


(Bl) 


1 - Ass(s) 

where ct; = 1 has been chosen as length unit and^® 

Lu{s) = L[°^+L[ps, (B2) 


LUs)=L\^J+L\1^s + L\^Js^ 


(B3) 


r (2) _ r (2) _ 

^Is — ^sl — 


1 l + g + 


3i7s 

l-''7i. 


12 4Tsiq^ -b77s/g 


(B12) 


This closes the determination of Gu{s) for given values 
of 77 s, g, and Then, by numerical inverse transform 
one can easily obtain gu{r). On the other hand, pure 
analytical expressions are also possible for the different 
layers l<r<l-|-g, l-|-g<r<l-|- 2 g, l-|- 2 g<r< 
1 -|- 3g, .... The trick consists in formally attaching a 
bookkeeping factor e to any exponential in G/;(s). Then, 
by expanding in powers of e we can write 


Gu{s) = Y, e-(i+"«)W„(s), (B13) 

n=0 


where we have made e = 1. From Eq. (B13) we get 


.j 0\-/ 

guir) = - Y - 1 - nq), (B14) 

n—0 

where "fnif) is the inverse Laplace transform of r„(s). 
The functions 7 n( 7 ’) can then be expressed in terms of 
the three roots of a cubic equation, analogously to the 
case of Eq. (6.2). Therefore, if we are only interested in 
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FIG. 11. Loci in the plane TsI vs rja where the PY approxi¬ 
mation predicts gii{r) = 0 at r = (1 -|- q)^ for q = 0.12 (upper 
curve) and q = qo (lower curve). The radial distribution func¬ 
tion gii{r) is not positive definite below each curve. 


the interval 1 < r < 1 -|- fcq, we just need to keep the first 
k terms in the sum of Eq. (B14). 

From a practical point of view, it is sufficient to deter¬ 
mine gii{r) in the interval 1 < r < 1 -|- 3g, in which case 
only 7o(r'), 7i('r), and 72 (f) are needed. Their analyti¬ 


cal expressions are easily obtained with a computational 
software program but are too lengthy to be reproduced 
here. In general, 71 (0) 7^ 0, what implies a jump discon¬ 
tinuity of gii (r) at r = 1 -b q, 

Sgii = gii{{l + q)~) - gii{{l + q)+) = 

1 + g 

^ (l + g)7?s [1 + q + {2 - q)T]s]'^ 

I2q{l - ?7 s )2 [(1 -b q)r]s + 4qTsi{l - ris)f 

(B15) 

Note that r = 1 -b g is the threshold distance beyond 
which no bridges are possible (see Fig. 2). This is clearly 
reflected by a strong decrease of gii(r) when going from 
r = (1 -b q)~ (bridges are possible) to r = (1 -b q)'^ (no 
bridging effect). This physical phenomenon can give rise, 
as an artifact of the PY approximation, to a negative 
value of gii(r) at r = (1 -b q)~^ if rjs is sufficiently large or 
Tsi is sufficiently small. This is illustrated in Fig. 11 for 
q = 0.12 and q = go- 

Once gii{r) is known, Eq. (3.2) gives the pair potential 
vii{r) in the PY approximation, as depicted in Fig. 5. 
The effective second virial coefficient can also be deter¬ 
mined analytically as follows: 


roc Q 

= -27r / dr [giiir) - 1] = 27r lini — [Gzi(s) - s“^] 

Jo 

8 -b (20 - 15g - 6g2 - g3)ry, -b 2(4 - 6g -b ?>q^ + q^Ws " 


/o 

TT 


12(1 + 277,)2 

2(1 -b g)?7s(l + q + 2gs- qVs) [6(1 + 2gs) + g^(l - ds)^ + 2g(l - ??s)(2 + Vs)] 
(1 - Vs) [(1 + q)Vs + ^qTsi{l - Vs)] 

(1 -b g)^7s(l + q + 2vs- qVs)‘^{2 + q + Ays - qVs) \ 

{l-Vs)[{i + q)Vs+4:qTsi{l-Vs)f J 


(B16) 


Appendix C: Derivation of the relationship between rjs, Vi, 
and in the SAO model 


where 




f dr^‘ ■ ■ ■ 


J dr^‘ 


(C3) 


In the semi-grand-canonical ensemble (zs,Ni, V, T) the 
average number of small particles can be obtained from denotes a canonical average over the bare solutes. Then, 


the associated thermodynamic potential B as 

d(3F 


{d'^s) Zs,Ni — Zs' 


dzs 


(Cl) 


taking into account that Q oc Zg, Eq. (Cl) reduces to 


Now, from Eq. (4.2) and the equality Uf^ = Uu + 11, we 


can write 


(Ns)z,M - - ■ 

Next, if g < goi ^ = ^0 + + ^2, so that 


{Ns)z^,Ni = —/dHo — — 






-pUn 


^g-/3n2) 


(C4) 


(C5) 


N, 




(C2) 


Note that the last term on the right-hand side can be 
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rewritten as 




drPvii{r)ges{r\rii,r]i^'>), (C6) 


where 




f dr<‘> ■ ■. f d4y'><- 


(C7) 


and we have taken into account that Ni{Ni — 1) ~ N'^ in 
the thermodynamic limit. Finally, applying Eqs. (4.7), 
(4.11), and (C6) in Eq. (C5), it is easy to obtain Eq. 

(4.20) . 

The first two terms on the right-hand side of Eq. 

(4.20) can also be obtained from the canonical ensem¬ 
ble {Ns, Ni,V,T). Up to the level of the second virial 
coefficient, the free energy F is 

^ = p,, In (psA^s) + Pi In {piAf) - Ps - Pi + ^-pf 

+2BsiPsPi + 0{p^), (C8) 

where Bsi = (1 + — 4^) ■ The solvent chem¬ 

ical potential is Ps = [d{F/V)/dps]p^, so that 


aPfd-s 

Zs = .3 = Ps [1 + 2Bsipi -b 0(p^)] , (C9) 


which is consistent with Eq. (4.20). 


Appendix D: SWAO model 

In the SWAO model Eq. (2.10) is replaced by 

{ - 1 , r < ash 

- 1, CTsi < r < (Tsj(l-b As/), (Dl) 
0 , r > CTs/(l-b As/), 

where tsi and UsiNsi are the depth and width, respec¬ 
tively, of the attractive well. One can define an effective 
stickiness parameter= 12 — l) As/, so that 

the SWAO model reduces to the SAO one in the double 
limit Cs/ —>■ 00, As/ —?► 0 at fixed Ts/. 

All the steps in Sec. IV up to Eq. (4.10) are still valid 
for the SWAO model. However, the condition for having 
= 0 if n > 3 is not CTs/ < cr/(l -b go)/2 (or q < qo) 
but fTs/(l + As/) < ai{l + qo)/2, what is equivalent to 
9(1 + As/) + As/ < qo- 

To simplify the expressions, in this appendix we take 
again cr/ = 1 as the length unit. Inserting Eq. (Dl) into 
Eq. (4.10), one obtains 

{ 00, r < 1, 

ip{r), 1 <-r < (1-b g) (1-b As/), (D2) 
0 , r > (1 -b g) (1 + As/), 


where the function ip{r) can again be decomposed into 
three terms (depletion-bpulling-bbridging), as given by 
Eq. (4.13), except that now 

ipp{r) = —[C(r-,l-bg)-bC(r+,(l-bg)(l + As/)) 

-2C{r,l + q)], (D3) 


iphir) = — [C{r,l + q) -bC(r, (1 -b g)(l -b As/)) 

-C(r_, 1 -b g) - C(r+, (1 -b g)(l -b As/))], (D4) 

where 

C{r,a) = ^{a — r)^{2a + r)Q{a — r) (D5) 

is the volume of a spherical cap of height a — r in a sphere 
of radius a and 


r± = r ± 


(1 + g)^ 

2r 


As/ 



(D6) 


The depletion term is still given by Eq. (4.14), i.e., 
tpd{r) = -(3/27rg3)C(r, 1 -b g). 

The ranges of the contributions ^d(^), 

'iphir) are 1 + g, (1 + q )(1 + As//2), and (1 + g)(l + A^/), 
respectively. It can be easily verified that in the sticky 
limit Asi —>■ 0 the potential of Eq. (D2) reduces to the one 
of Eq. (4.12). One can also verify that the jump discon¬ 
tinuity at r = 2cFsi= 1 + g of the SAO model disappears 
in the SWAO one, which is everywhere continuous. 
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